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1. Introduction 

Matrix product states (MPS) have attracted considerable interest in the interdisciplinary 
field of condensed matter physics and quantum information science [HIS]- The archetype 
of MPS can be traced back to the seminal work of Affleck, Kennedy, Lieb, and 
Tasaki [31 HJ |5], which is aimed at understanding the ground state properties of 
Haldane gap systems O [7]. It was then formalized in more generalized and abstract 
ways [HllHllin]. After those works, the density matrix renormalization group (DMRG) 
method, a powerful numerical method to study one- dimensional strongly correlated 
systems, was introduced by White [HI [I2] and its connection to MPS formulation was 
discussed by Rommer and Ostlund [13]. A common feature of mathematically rigorous 
and numerical approaches is that both have failed to describe quantum critical ground 
states which exhibit quasi-long range order. In the context of quantum information 
theory, Vidal et al, have attempted to characterize the quantitative difference between 
non-critical and critical ground states in one dimension in terms of entanglement 
entropy [14]. However, a natural question to ask is whether or not there are quantum 
critical ground states expressed in a form of MPS. Surprisingly, the answer is yes. In 
Ref. [15] and the subsequent work [16], Alcaraz and Lazo have actually shown that 
the eigenstates of the spin-i Heisenberg chain can be expressed as MPS. It has been 
well known that this model is gapless and the quasiparticles, so-called spinous, have a 
linear dispersion relation since the pioneering work of Bethe in 1931 [17]. The method 
to solve this model used by Bethe is called coordinate Bethe ansatz and there are 
several variants (algebraic, functional, etc.,) of it. Alcaraz and Lazo proposed another 
alternative formulation of the Bethe ansatz in terms of MPS, which they call matrix 
product ansatz (MPA). In this formulation, one can obtain the spectrum conditions, i.e., 
Bethe equations, imposing algebraic relations between matrices consisting of MPS. The 
physical meaning of those matrices is interpreted as the generators of the Zamolodchikov- 
Faddeev algebra in (1+1) dimensional field theories [18]. In Ref. [19], Alcaraz and Lazo 
applied MPA to other integrable models such as the Hubbard model [20], fermionic 
supersymmetric t-J model [21], and Fateev-Zamolodchikov model [22], and obtained 
the correct Bethe equations. 

In this paper, we show that MPA is essentially equivalent to the algebraic Bethe 
ansatz (ABA) in the XXX and XXZ spin-^ Heisenberg chains. The ABA is an elegant 
method for solving the eigenvalue problem of quantum integrable models developed in 
the late 70s [231 1211 [25] • So far, the relation between the MPS and ABA has been 
discussed in a completely different context, i.e., stochastic Markovian models in one 
dimension [26] . One of the simplest examples of the models is the asymmetric simple 
exclusion process (ASEP), which plays an important role in non-equilibrium statistical 
mechanics. This model was first exactly solved not using the Bethe ansatz [27] while 
the relation to the integrable quantum spin chains was then clarified by Alcaraz et al ., 
[28] . It was first discovered in Ref. [26] that the MPA for ASEP can be derived from 
the ABA. The authors have also constructed the explicit finite dimensional matrices 
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for the MPA. The key ingredient is the change of basis in J^, the tensor product of 
auxihary spaces. Since the MPS are defined as traces over M' ^ they are invariant under 
the change and hence one can take an appropriate basis in which the matrices have very 
simple expressions. Along the same lines as their approach, we derive the MPA for the 
XXX and XXZ spin-^ Heisenberg chains. The explicit expressions for the matrices are 
also obtained. 

The organization of this work is as follows. In Section 2, we review the ABA 
method for solving the eigenvalue problem of the Heisenberg chain. In section 3, we 
derive the MPS representations for the Bethe eigenstates from the ABA by preparing 
J^, the tensor product of the auxiliary spaces. In section 4, the explicit expressions for 
the matrices appearing in the MPA are obtained by solving the recursion relations for 
the matrices. In section 5, the relation between the MPA and the six-vertex model with 
domain wall boundary conditions [22] is discussed. Conclusions and future perspectives 
are given in the last section. In Appendices, we provide graphical representations for 
the matrices in the main text and discuss the mapping from the six-vertex model to a 
five-vertex model. 



2. Algebraic Bethe ansatz for the XXZ spin-| Heisenberg chain 

The XXZ spin-i Heisenberg model with the periodic boundary condition is described 
by the following Hamiltonian: 

L 

Hxxz = + ^'^'+1 + - 1)} ' (1) 

i=\ 

where L denotes the total number of sites and af [a = x, y, z) are the Pauli matrices 
defined on the ith site. Here A is the anisotropy parameter and the particular cases 
A = 1 and A = correspond to the XXX and XY chains, respectively. 

The eigenstates of this model can be constructed using the ABA [23l [2ll [25l [30] . 
We shall briefly review this construction. In the ABA, the central object is the quantum 
i?-matrix which is the solution of the Yang-Baxter equation. For the XXZ model, the 
i?-matrix acting on (g) is given by 

/ sinh(A + r/) \ 

sinh A sinh rj 

sinh?7 sinh A 

y sinh(A + rj) J 

where A is the spectral parameter and the relation between rj and the anisotropy A 
is given by A = cosh?]. Next, we introduce the quantum L-operator represented by 
a matrix acting on the tensor product of two-dimensional vector spaces Vj ® Vi. The 
auxiliary space Vj introduced here is spanned by two orthonormal states labelled as |<— ) 
and 1—7-) while the physical Hilbert space at the ith site Vi is spanned by |t) and ||). 
For the XXZ model, the L-operator Cji{X) is defined by the i?-matrix as 

£,,(A) = R,, (a - I) , (3) 



RiX) ^ 



sinh(A -|- 77) 



(2) 
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where Cji acts on Vj ®Vi. Note that the operator Cji{X) acts trivially on all the sites 
other than Vi. More explicitly, the L-operator is written as 

/ 1 \ 

b{X) c(A) 

c(A) b{\) 

\ 1 / 



(4) 



where 



^^^^ - sinh (A + ' '^^^ - sinh (A + 2) 



for the XXZ chain and 

"W-^. ^(A)~ (6) 

for the XXX chain. The latter can be obtained by taking the scaling limit of the former, 
i.e., Axxz = eAxxx and r] = ie with the limit of e — )■ 0. It is useful to note the relation 
between 6(A), c(A) and the physical quasi- momentum k. For both the XXZ and XXX 
cases, 6(A) = z with z = e~^^ [15]. By a straightforward calculation, one can verify the 
following relation: 

c(A)2 = 1 - 2A2 + z\ (7) 

In particular, c(A) = 1 — z when A = 1 corresponding to the XXX chain. The following 
identity can then be shown as a direct consequence of the Yang-Baxter equation: 

Roi{\ - /i)£oi(A)£ii(/i) = £ii(/i)£oj(A)i2oi(A - /i). (8) 

Here, i?oi acts nontrivially on Vq® ^1 and trivially on Vi. Note that Cu acts nontrivially 
on Vi ® ^ and trivially on Vq and the other spatial sites. 

The monodromy matrix is then constructed as the following ordered matrix 
product: 

ro(A) = £oi(A)/:o2(A)...£oL(A). (9) 
In the basis of Vq, the monodromy matrix can be represented as a 2 x 2 matrix: 

where matrix elements A(A), B(A), C(A), and D(A) are themselves operators acting 
on the total Hilbert space = ®^iVi. Using the relation Eq. ([H]), one can show the 
following relation for the monodromy matrix. 

i?oi(A - MiXmifi) = ri{fi)%{X)Roi{X - /i). (11) 

Graphical representations for both Eq. ([8]) and f|TT]) are shown in Appendix B.2. The 
commutation relations among A, B, C, and D can be obtained from this relation. 
Taking the trace of the monodromy matrix over Vq, one obtains a one-parameter family 
of transfer matrices acting on J^: 

T(A) = Trviro(A) = A(A) + D(A). (12) 
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The Hamiltonian Eq. ([T]) can then be obtained from T(A) by the trace identity: 

d 

Hxxz = 2 sinhr/ — log T(A) + const. (13) 

OA X=r]/2 

Since the Hamihonian commutes with the monodromy matrix, one can construct a 
simultaneous eigenstate of both Hxxz and T. The eigenstate of this operator is 
constructed by B(A)'s as 

|Ai,A2...A„) = B(A„)...B(Ai)<l>o, (14) 

where $o denotes the reference ferromagnetic state, i.e., $o = |lT) = I t^t ••• t) and n 
denotes the number of down spins. We call the above state a Bethe state. Since one 
can show that B(Aj)'s commute each other from the relation ( |TT1) . this state is invariant 
under permutations of Aj's. The spectrum conditions, i.e., the Bethe equations, are then 
obtained using the commutation relations between A, D, and B [231 [2^ 125] . For the 
XXZ model, those equations are given by 

(^^l!ra)'n "ta '^^'^1 =^- ^=i.2.-."- (15) 

Vsmh(Aj + j) J smh(Aj ~ Xk - V) 

3. Derivation of the matrix product state representation from the algebraic 
Bethe ansatz 

In the previous section, we outlined the construction of the eigenstates of Hxxz using 
the ABA. In this section, we derive the matrix product state representations for the 
eigenstates from the ABA. Let |-(— ) and |— )■) be two orthonormal states spanning V . 
Then the Bethe state is expressed as 

|Ai,A2...A„) = B(A„)...B(A2)B(Ai)$o 

= (^1 T{\n) K) ... (^1 r(A2) K) (^1 r(Ai) k) $0 



Try®n(g„r(A„) ® ... ® r(A2) ® r(Ai)$ 







Tr y®n Qr 



$0 (16) 



with Qn = K, ... ^ ... ^1 = I ^)(^ I and ^(Ai, A„) = £i(A„) ® ■ ■ ■ (g) 

£j(Ai) which acts on the vector space Vi only in J^. Here, we have omitted the indices 
for the auxiliary spaces. It would be helpful to note that the following identity for tensor 
products holds: (A.B) ® {CD) = {A ® C).{B ® D). 
We now introduce two matrices Dn and C.„ via 

/:,(Ai, A„)| t) = A.(Ai, A„)| t) + C„(Ai, A„)| i) (17) 

or equivalently, ^^(Ai, A„) = (t |£i(Ai, A„)| t) and C„(Ai, A„) = (| 
|£j(Ai, A„)| t). It should be noted that Dn, Cn, and Q„ are 2" x 2" matrices acting 
on V^®" with scalar elements. The first terms are, for example, given by 
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Henceforth, we shall denote the n-fold tensor product of the auxiliary spaces as 
M' = <S>^=i ^n-j+i- Using Dn, Cn, and Qn, the Bethe state can be written as 

L 

|Ai, A2, A„) = Tr Q„ lliD^Xi, A„)| t). + C„(Ai, A„)| , (19) 

where | t)i and I i)i denote the up and down spin states at the zth site, respectively. 
This form can be regarded as the matrix product state in the usual sense except for the 
boundary matrix Qn- 

Let us consider the recursion relation between Dn+i, C„+i and D„, C„. From the 
fact £i(Ai, A„, A„+i) = Ci{\n+i)®Ci{\i, A„), one can derive the following relations: 

L'„+i(Ai...A„,A„+i) = I I ^ ) ® D„,(Ai...A„) + f ° ° I ® C„(Ai...A„,) 

U 0[An+l) I \ C[An+l) U 



c(An+l)C„(Ai...A„) b{Xn+l)Dn{\l...Xn) 



(20) 



and 



C„+i(Ai...A„,A„+i)= I Q '^Y'M ( ^^V'^ 1 \^Cn{Xi:.\n) 



(21) 



b{Xn+l)Cn{Xi...\n) c( A„+i )-D„, ( Ai . . . A„ 

a(Ai...AO 
Here we have used the following identities: 

.(tiA(A„«)it).=i ; ,,,°^^)), .(tiA(A„«)ii).=(^„°^^) I 

,(i |£.(A„«)I t), = I ° "'^S*'* ) , |£.(A„«)I i)i = f ''^S"'* I 

The above matrices are understood to act on the auxiliary space Vi (see Appendix B.l 
for a more detailed explanation). By definition, it is obvious that (5n(Ai...A.„, A.„+i) 
satisfies the following recursion relation: 

Q.+i(Ai...A„,A.+i)= 0) ^^^^ 

and hence the explicit form of Q„(Ai, A„) is given by 

QniXl,...,\n) = 0(^[ (23) 

Therefore, one can construct a matrix product representation of the Bethe ansatz state 
using the relations fl20tl22p recursively. The interesting point here is that the dimension 
of the matrix is finite even in the thermodynamic limit if the number of down spins is 
finite. On the other hand, if we consider the case of the fixed magnetization, i.e., the 
ratio n/L is fixed, the dimension of the matrix becomes infinity in the limit of L — )■ 00. 
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4. Change of the basis in 

In the previous section, we have derived the matrix product state representation for the 
Bethe states from the ABA. At that point, however, the matrices are defined by the 
recursion relations and the connection to another representation proposed by Alcaraz 
and Lazo [T5l [16] is not clear. In this section, we clarify the direct relation between 
them by solving the recursion relations. We now try to rewrite the matrices and C„ 
as diagonally as possible. For simplicity, we omit here the indeterminants of matrices, 
i.e., Ai,...,A„. Furthermore, we will replace h{\j) with Zj = e"*'^^. Suppose that -D„ is 
diagonalized by the invertible matrix Fn as F~^DnFn = X'n- In the new basis in J^, Cn 
and Qn are transformed into C„ = F~^CnFn and Q„ = F~^QnFn, respectively. From 
the cyclic property of the trace, it is obvious that the Bethe state can be written by new 
matrices as 

L 

Qnl[iVr,.\t)^ + Cn\i)i 



(24) 



|Ai, A2, A„) — Tt 
We now decompose C„ as a sum of n matrices, 

n 

Cn = 5^C«, (25) 



i=l 



and suppose that P„ and Cn^ satisfy the following algebraic relations: 

C«I)„ = z,VnC^^ (26) 

C^^Cl^^ = S.,{z.,z,)Ci;^Clp (27) 

C»C« = 0. (28) 

Then we look for an appropriate Sij{zi, Zj) by observation. Let us now consider the case 
of n = 2. In this case, D2 and C2 are given by 

D = I ^'^^'^ ^ ] + ( ° M (2^) 

' ' Z2D^i\^) )^[ c(A2)Ci(Ai) J ' ^^^^ 

, ci\2)D^i\^) \ ( z^dX,) A 
- ' Ci(AO j' 

respectively, where Di{\i) and Ci(Ai) are defined in Eq. (ITS]) . By an explicit calculation, 
one can confirm that the matrix D2 is diagonalized as 

V2 = F^'D2F2= \ ^'^^'^ ° , ), (31) 



where 

with J^i = ^^z^C'i(Ai). We note that the inverse of F2 is given by 



(32) 




F.-' = V . (33) 
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In the new basis defined by F2, the matrix C2 becomes 

-.,c(A,)^, + C,(AO 

We now divide the above matrix into two matrices as C2 = C^^ + C^^ with 



(34) 



f Z2C,{X^) + ciX2)Ti \ (2)_ / c{X2)D,{X,) \ 

' ~[ -c(A2)^iJ-i + Ci(Ai) ^ " 1^0 j-^''^ 

One can confirm that Cg^^ and Cg^'' satisfy the algebraic relations (!26|) and (128|) . Then 
we compare the commutation relation between C^^ and C^'' with Eq. ( l27j) and find that 
5*12 (-21, Z2) should be given by 

Q f \ ^1 Z1Z2 + I-2AZ2 , , 

From this observation, we take the coefficient Sij{zi, zj) to be 

S- (z- Z-) - ■ '''' + ^ ~ f37) 
z,z, + l-2Az, ^^^^ 

for general i and j. We shall prove the algebraic relations ( I26ll28l) with this Sij{zi, Zj) by 
induction on n. For n = 1, the matrices Vi = Di, Ci = Ci trivially satisfy the relations 
( l26ll28l) . Now, we suppose that we have already found the matrix F„ which diagonalizes 
Dn as T>n = F~^DnFn and have found a decomposition C„ = J2^=i "* ^^'^^ that the 
relations f l26H28|) are satisfied between and Cn'^'s. Then we show that it is possible 
to construct a matrix Fn+i which diagonalizes Dn+i and can find a decomposition of 
Cn+i- First, we take the matrix -F„+i to be of the following form: 

Fn^^=\ Jl I 1 (38) 



F -F F 



or equivalently 



We then obtain the following recursion relations using Eqs. f l20H22|) : 



^ ~J'rJ^n + c(An+i)C„ + Zn-\-lT>nJ' n ^n+l^^n ^ 

Q _ I ^n.+lC„ + c(Xn+l)'DnJ' n c(An+i)!Dn . ^^^-^ 

^n-F n Z,i2-\-lF nCn c(Aj^-|_i)F" x\T^nJ~ ^ Cn c(^Xn-\-l)F riT-^n 



Qn+l = \ ^ , , (42) 




where Fn is an unknown matrix to be determined. Since one can take Fi to be the 2x2 
identity matrix, Qi = Qi holds and hence Qn = Qn for ^n. In other words, the change 
of basis preserves the domain wall boundary condition (DWBC) in J^. This property 
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plays a crucial role when we interpret our results in terms of the six-vertex model as 
discussed in Sec. 5. From Eq. f HOjl . one can see that P„+i is a diagonal matrix if the 
matrix Tn satisfies —J^rJ^n + c{\n+i)Cn + Zn+iDnJ^n = 0. Using Eq. ( 126|) . such J^n can 
be constructed as 

J-. = c(A„+ip-iJ]--^. (43) 



Next, we shall show that the (2,l)-component of Cn+i is also zero if Sij{zi,Zj) is given 
by Eq. fl37|) . By a direct calculation, one obtains 

n n 

= ciK+lWn'Y.Y.'^^fnC^n (44) 
i=l j=l 



with 



1 1 1 c(A„ 



dij Zj^ -^71+ 1 



\2 



~ Zn+1 Zi — Zn+1 {Zi — Zn+l){Zj — Zn+l) 

_i 1 1 1 - 2A2;„+i + 

\ ^n+l 



Zj — Zn+1 Zi — Zn+1 {Zi — Zn+l){Zj — Zn+l) 

Zn+1 ZiZj + 1 — l!S.Zi 



(45) 



i^i ~ Zn+l){Zj — Zn+l) 

where we have used the relation Eq. ((71). Using the relations ( 1271) and ( |28|) . Eq. (14^ is 
zero if dij satisfies dij + [Sij{zi, Zj)]~^dji = 0. Although this relation is highly nontrivial, 
one can confirm that it holds if Sij{zi, Zj) is given by Eq. fl37|) . The recursion relations 
Eqs. ( l40l) and ( 14T1) can now be written as 

Vn+i= { I ° I (46) 

Zn+l J 

^^^ ZiZ„+i+l-2Az„+i 



En ZiZ„+i+l-zn.z„+i /^(i) / \ 

i=l 1:^1:7, C[An+l)^^n 

»+! + 

^n + 1 

n+1 p{i) 



C — ^«=i z^-z„+i (AJ\ 

'-'"+1 I n -V Y^" z,z„+i+l-2/^Zi \-\^') 

U -Zn+l 2^i=i (-n 



From them we deduce the decomposition of Cn+i = Yl^=i Q+i with 

+ 1 - 2AZn+l)Ci'^ 
""^^ -Zn+l{ZiZn+l + l -2AZi)C, 



and 



^(n+l) _ f c{Xn+l)Vn \ , . 

Cnn-i - I Q Q j- (49) 

From the above expressions, it is shown that Vn+i and {I < i < n + 1) also satisfy 

the algebraic relations Eqs. ( I26|) -(l28l). Therefore, we have proved that there exists a 
decomposition of C„ with algebraic relations ( !26ll28l) at any n. 
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We now write down the explicit expressions for D.„ and Cn^ using the recursion 
relations ( I46ll49l) . From Eq. ( 146|) . we first obtain 





-i+i 

Then, from Eq. ( H9l) . we obtain 



(50) 



We finally find the explicit form of Ci* ^ with 1 < i < n. The relation f l48p can be recast 



as 



where 



pit) _ ^i^n+l + 1 — 2AZi I Si^n+l{Zi, Zn+l) \ ^ pii) / r:o\ 



X 1 - 2A^- • ^^^^ 
Therefore, we obtain the explicit expression for as 

c« = /(.„...,.„)c(A.)(^| ^g,^o 1^ ^ ^1 

with 



/(........„)= n (55) 

l=i+l ^ Zi Zi J 

for 2 = 1, n — 1. We note /(^n) = 1- The action of the matrix is almost diagonal 
in ^ except in the single auxiliary space Vi. Therefore, it is now quite trivial to confirm 
the algebraic relations fl26l l28|) using Eqs. fl50|) . fl5T]) and (154]) . 

Let us now clarify the relations between the ABA and the MPA explicitly. Keeping 
the notations as close as possible to those in [T5| [T6] , the Bethe eigenstate in the MPA 
is written in the form: 

iV'n) = 5^ Tr {El^-^A^El^-'^^-^Ar, ■ ■ ■ X2, (56) 

{X1,X2,...,X„} 

where |a;i,X2, ■■■iXn) (1 < xi < X2 < ... < Xn L) denote the configurations with down 
spins at (xi, X2, x^) and the subscript n indicates the total number of down spins. 
Note that in the original papers by Alcaraz and Lazo [T5l [16], (xi, X2, x„) are the 
locations of up spins. The matrix An is decomposed by n matrices as 

n 

= ^ ^ Ak^^nEyi, (57) 
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where the matrices A^.^n obey the commutation relations 

A-ki,nEn = ZiEnAk-^n, (58) 
Aki,nAkj^n — Sij{Zi, Zj)Akj^nAki^n: (^9) 
EMn = e-'^'rinEn, (60) 

with P = Yl^=i^i- The above relations assure that \ipn) is the eigenstate of the 
Heisenberg Hamiltonian Hxxz- From Eq. ( 159|) . one can derive A^.^^ = 0. The relations 
( I58ll59|) together with Al^^ = look very similar to Eqs. fl26tl28p . We now try to find 
a one-to-one correspondence between the matrices in Eqs. ( I26ll28l) and those in Eqs. 
( I58ll59l) . To reproduce the correct commutation relations, the following relations are 
required: 

En = aVn, = /3Qn = /3Qn, and Ak,,nEn = 1^C^\ (61) 

where a, /3, and 7^ {i = 1, ...,n) can be the arbitrary numbers. Here we have used the 
fact that ZjSij{zi, Zj) = ZiSij{zi, zj). For simphcity, let us fix a = /3 = 1. Then, from 
the above correspondence, we obtain Ak^^n = liCn'T^n^ and find 

A,,,, = i^nz,, .....„)(§( ° 

^* 1=1 \ 

(63) 

which are indeed equivalent to the matrices found in [16] if we change the ordering 
of quasi-momenta from {ki, k2, kn) to (/c„, /c2, /ci) and take an appropriate set of 
7j's. In this way, we have derived the matrices appearing in the MPA using the ABA. 
Therefore, the matrix product Bethe ansatz is equivalent to the algebraic Bethe ansatz. 
We note here that if we take the XY limit (A = 0), S{zi, zj) = —1 and hence Eq. (163!) 
can be regarded as the Jordan- Wigner transformation in the auxiliary space J^. 

We remark on the coefficient 7^ in the proportional relation between A^^^nEn and 
Cn\ In the original work by Alcaraz and Lazo, An is defined by Yl^=i ^ki,nEn = 
Yl^=i'yi^n\ while Cn is defined by Eq. ( I25l) . Therefore, one may think that there 
is a constraint on the proportionality coefficients 7j's. However, they can be arbitrary. 
Let us explain the reason for it. Since A^. „'s are nilpotent {Al, ^ = 0), only the following 
products appear in Eq. fl56l) : 

where a are permutations of (1,2, ...,n). Therefore, if we take an arbitrary set of 7j's, 
the matrix product Eq. (164|) for any a has the same prefactor and hence the state is 
uniquely determined apart from an overall factor. 

We also remark on the relation between Eqs. fl24|) and fl56|) . Using the 
correspondence (|6T1) . we can rewrite Eq. as 

L 

nnl[{En\t)i + An\l)^) , (65) 



|Ai, A2, A„) — Tr 
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where we have set a = /3 = 7^ = 1. It is easy to show |Ai, A2, A„) = IV'n), where lifj^) 
is defined in Eq. f l56|) . One may think that in the expansion of the r.h.s. of Eq. fl65|) 
there are traces of the following matrix products 

where {11,12, im) are subsets of (1, 2, n). Such terms do not appear in Eq. ( 156|) . In 
fact, they are forbidden since Tr ^-[fi„M„] is nonzero only when the matrix M„ is written 
as Mn = + ■ ■ ■) where denotes matrix transpose. Let us explain it in more detail. 
From Eq. ( 163|) . it is obvious that only a single arrow which is one of the orthonormal 
vectors in Vi is fiipped by the action of A^^^n- On the other hand, Tr^[r2„M„] can be 
recast as {-^ |M„| =^). This matrix element is nonzero only when all the arrows are 
fiipped by the action of M„. This proves that the trace of (l66l) is zero when m < n. 



5. Relation to the six- vertex model with domain wall boundary conditions 

In this section, we shall clarify the relation between the matrix product Bethe ansatz 
and the six-vertex model with domain wall boundary conditions. The six vertex model 
is a two-dimensional statistical mechanics model in which the Boltzmann weights are 
assigned to the six different configurations of arrows around a vertex. If the Boltzmann 
weights satisfy the Yang-Baxter relation, the model is exactly solvable by the Bethe 
ansatz. The partition function of this model on a n x L rectangle is defined by 

1 ^ 



3f= II sinh ( A, + ^ ) z with z=J2 U [A(A,)](;:;:, (er) 

config v={i,j) 

where the summation is taken over all the possible configurations satisfying the ice rule 
and V = denotes the vertex which is the intersection of the zth vertical and jth 
horizontal lines. In the definition of Z, the product is taken over all the vertices. The 
Boltzmann weights are related to the six nonzero matrix elements of the L-operator and 
Ci{Xj) is assigned to the vertex v. The indices fi and u correspond to arrows on the 
horizontal edges while a and p correspond to spins on the vertical edges (see Fig. [I](a)). 
The model with periodic boundary conditions was first solved by Lieb [21] • The model 
with domain wall boundary conditions (see Fig. [1] (b)), which is relevant to our study, 
was first discovered by Korepin in the context of the norm of the Bethe state [22]. For 
a finite (n x n) lattice, the partition function apart from an overall factor is given by a 
component of the Bethe state as 

^r^^ = (^|Ai,...,A„), (68) 

where (JJ. | = i{i \- In Ref • [29] , the recursion relation for Z^^^^ was derived, 

and then later it was solved by a determinant formula [32] • Note that in the original 
papers, an inhomogeneous generalization of the six- vertex model was studied. Not only 
the component of the state | but also the other components of the Bethe state can be 
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(a) 



a 



(b) 



V \' V 



(C) 



Figure 1. (a.) L operator associated with the vertex which is the intersection of the 
ith vertical hne and jth horizontal line, (b) Domain wall boundary conditions, (c) 
Other boundary conditions. Each ai is or 4- 



interpreted as partition functions of the six-vertex model. The coefficient of the state 

(a!^^l,...,af^l|Ai,...,A„) (69) 

corresponds to the partition function with the boundary conditions shown in Fig. [T] (c) 
(A more detailed discussion is provided in Appendix B. 1). Here, cxj {i = 1,2, 
denotes the spin state, i.e., t or I. In this sense, to obtain the Bethe state is equivalent 
to obtaining the partition function of the six-vertex model with various boundary 
conditions. It appears to be a formidable combinatorial task. However, if we construct 
a vertex model corresponding to the new basis in J^, this complexity is greatly reduced 
as shown in Appendix B.4 and 5. It corresponds exactly to the similarity transformation 
from the triangular matrices Dn and C„ to the matrices Vn and Cn^ whose actions are 
diagonal and almost diagonal except for the single space Vi, respectively. It is important 
to stress here that this change of basis does not alter Qn, i.e., the domain wall boundary 
condition in J^. In the corresponding vertex model, the Boltzmann weight of one of 
six vertices becomes zero. We henceforth call this model a five-vertex model. Note that 
since the partition function itself is the same after the transformation, the statistical 
model is unchanged. Compared with the six-vertex model, however, the combinatorial 
complexity has been already resolved in the five-vertex model and hence the possible 
number of configurations involved in the calculation of Z is greatly reduced. 

A similar reduction has been known as F-matrices in the context of the Drinfel'd 
twist and triangular Hopf algebras [331 El] which have been used to calculate the 
correlation functions in the Heisenberg chain [351 EHl EZl |38]. The F-matrix is named 
after a factorization of the i?-matrix and defines a transformation from A, B, C, and D 
to the new operators. By this transformation, A and D are transformed into diagonal 
matrices. On the other hand, B and C are transformed into the matrices whose actions 
are almost diagonal in J^. By expanding the physical spin operators in terms of those 
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diagonal and almost diagonal operators, it is able to handle calculations of various 
correlation functions. Compared with this approach, the matrix F„ we found is a 
counterpart of the F-matrix in . Note that the original F-matrix is defined in the 
physical space M' while our F„ acts on M' . It is surprising that the connection between 
the MPA and the ABA we found turns out to be related to the F-matrices. It would 
be interesting to bridge over these two approaches more concretely. 

6. Conclusion 

In conclusion, we have derived the matrix product representation of the Bethe ansatz 
state for the XXX and XXZ spin-| Heisenberg chains from the algebraic Bethe 
ansatz. We have also shown that the finite dimensional representations for the matrices 
appearing in the matrix product ansatz proposed by Alcaraz and Lazo are equivalent 
to those obtained from the algebraic Bethe ansatz by use of the nontrivial change of 
basis in M' ^ which is related to the F-matrices. In the new basis, the matrices have 
a quite simple structure and the algebraic relations between matrices can be shown 
very easily. The relation between the MPA and the six-vertex model with domain 
wall boundary conditions has also been discussed. It would be of great interest to 
apply the obtained explicit matrices to calculations of static and dynamical correlation 
functions [39l|3ll[36l[371|38lll0lllllll2] and entanglement properties [l3l El iS] 
in the Heisenberg spin chains. The first step should be a calculation of the norm 
of the Bethe eigenstate, which has been obtained and is related to the determinant 
expression for the partition function of the six- vertex model with domain wall boundary 
conditions [291 |32l |16]. It would also be interesting to find a relation between the matrix 
product ansatz and the hidden Grassmann structure in the XXZ model discussed in 
the context of mathematically rigorous approaches |171 HHl HH]. Another interesting 
direction is a systematic construction of the matrix product state representation of the 
Bethe states for other integrable models, especially for correlated electron systems such 
as the Hubbard model [201 [501 [511 [52] . Our discussion in Section 2 indicates that it can 
be obtained if the model is exactly solvable by the algebraic Bethe ansatz. However, it 
is highly nontrivial whether or not we can obtain explicit expressions for the matrices 
since we need to find a new basis in which the matrices are almost diagonal. 
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Appendix A. index of notation 

Here is a list of tlie notations wliicli are used in tlie main text. 

• Linear spaces 

— V: a. pliysical Hilbert space spanned by |t) and ||) 

— V : an auxiliary space spanned by |— )■) and |^) 

— Jif: V^^ 

— jr : V^"" 

• States 

— 1^): |t,t...t) (=<^>o) 

— \^)- 1^,^ ... ^) 

— K): K,^ ... ^) 

• Parameters 

— b (or z) c : elements of L-operator 

— A, 77: parameters of Hamiltonian 

— A: a spectral parameter 

• Operators acting on the physical space {^) 

— A, B, C, D : elements of a monodromy matrix used in ABA 

— H : Hamiltonian matrix 

— T : a transfer matrix 

• Operators acting on the auxiliary space {J^) 

— C,D: matrices used in MPS 

— C, V: matrices used in the new basis after transformation 

— F : a. similarity transformation 

— : a matrix defined by Eq. (H3l) 

— R : i?-matrix 

— Q, Q : matrices for DWBC. (They turned out to be the same). 

— A, E: matrices used in MPA by Alcaraz and Lazo [151 HB] 

• Other operators 

— C : L-operator 

— T : a monodromy matrix 

Appendix B. Graphical representation 

Appendix B.l. definitions of graphs 

The L-operator Cji{Xj) given by Eq. (j4]) is a 4 x 4 matrix acting on a vector space Vj^Vi. 
The space Vj is spanned by | and | while Vi is by | t)i and | We hereafter 
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denote it as £j(Aj). The order of the basis in Vj<^Vi is fixed as |^t) > l^i) j l~^t) > 
In a graphical representation, we denote a matrix element as 



II — 



where /i and z/ are ^ or — )■ while a and p are t or |. The space Vi (or Vj) is denoted by 
a black thick line (a blue bold dashed line). There are 2^ = 16 possible configurations of 
vertices. However, 6 of them are nonzero as we have seen in Eq. which corresponds 
to the six-vertex model. The matrix Ci{Xj) is graphically represented as 



( 1 








o\ 







c(A,) 








c(A,) 


6(A,) 




























\ 















v 



(B.l) 



The correspondence between the vertices and the "Boltzmann weights" is summarized 
as follows: 



1, 



c. 



Note that the weights b and c depend on j and can be complex. We should stress here 
that the above correspondence is not unique but depends on an order of basis. We 
now introduce a graphical representation for the L-operator. The L-operator acting on 
Vj ® Vi is drawn as a vertex without arrows: 



A(A,) 



Next we define four kinds of operators acting on Vi from the L-operator. They are 
defined by decomposing the 4x4 matrix in Eq. (IB. II) into 2x2 sub-matrices: 



/j 



\ 



Here, the index \V~\ denotes the space in which the L-operator is represented as a 2 x 2 
matrix. In addition, we can define the operators acting on Vj by a decomposition of the 
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4x4 matrix: 



/ 1 








o\ 





6(A,) 


c(A,) 








c(A,) 






























V 



Ai 



Here, the index [V] denotes the space in which the 4x4 matrix is represented as a 2 x 2 
matrix. Although the 4x4 matrix appearing in the above equation is exactly the same 
as £j(Aj), i.e., the matrix in Eq. (IB.ip . the order of basis has been changed and hence 
the graphical representations for the matrix elements are different from Eq. (IB.ip . The 
reason why the 4x4 matrices themselves are the same is related to the fact that the 
Boltzmann weights of the six-vertex model are invariant under the simultaneous reversal 
of all arrows when there is no external field. Two operators out of the four operators 
acting on Vj are defined in Eq. ( 1T8|) . and graphs of them are given by 



ii 



Ai 



Note that the absence of arrows on the horizontal (dashed) line indicates that the above 
operators act on the space Vj. 

We are now ready to introduce a graphical representation of the Bethe state. As 
was seen in the main text, the Bethe state can be defined through the monodromy 
matrix T(Aj) or £i(Ai, . . . , A„). These operators are drawn as 



r(A,) = (g)/;.(A,) = n 



(B.2) 



/j 



n-1 



£j(Ai, . . . , A„) 



1=0 



n-1 

n 



\ 





j 








'i 


■--V 


j 








ii / 



2 
1 



(B.3) 



Here, each bond connecting two vertices represents summation over two possible arrows. 
The monodromy matrix is represented as a horinzontal object while £i(Ai, A„) is 
a vertical one. The connection between nearest vertices can be easily interpreted as 
the multiplication of 2 x 2 matrices. A creation operator for the Bethe state (B(Aj)) 
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corresponds to the (l,2)-component of T{Xj) which acts on 
represented by 



Y®L_ Using B{Xj 



B(A,)= {^\n\,)\^) = L- 










1 


2 





the Bethe state is drawn as 



. . . , A„) = B(A„) ■ ■ ■ B(Ai 



On the other hand, operators D„(Ai, . . . , A„) = (t| £j(Ai, . . . , A„) |t) and 
C„(Ai, . . . , A„) = {l\ Ci{\i, . . . , \n) It) on ^ = are represented as 




-Dn(Ai, . . . , A^ 



Cn(Ai, . . . , An) 



We shall draw a figure for each component of the Bethe state. From the matrix product 
representation given in Eq. ( lT6l) . each component is written as 

L 

Qnl[{(Ti\Ci{Xi,...,Xn)n)i 



/JVl] JVl] ^ X \ _ _ 



1=1 



n[(^.|A(Ai,...,A„)|t)JK). (B.4) 

1^), and 



i=l 



To draw a graph of Eq. flB.4p . we further introduce patches for ( 
((Til Ci{Xi, A„) It), in the following way: 



2 
1 



((Ti|A(Ai,...,A„)|t) 



ii 



It is now obvious that the graph for the component of the Bethe state |Ai, . . . , A„) is 
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represented as 



[a^ , . . . , Cr^ . . . , An) 



il i2 i AL 

The DWBC in corresponds to the case where VcTi =1. The partition function of the 
six-vertex model with this boundary condition is graphically given by 




7DWBC 



(^|Ai,...,A„) 




(B.5) 



which is the coefficient of the all down (full-filled) state. In particular, if L 
shown that the partition function is expressed as 

n-l 
1=0 



n, it is 

(B.6) 
(B.7) 



Appendix B.2. Yang-Baxter equation 

To illustrate the Yang-Baxter equation (YBE) in Eq. ([8]), we introduce a graphical 
representation for the i?- matrix which is defined by Eq. ([2]): 



J I 



(B.8) 



Here, the parameter Aj (or A^) is specified by the space Vi {Vj) for Rji{\j — A^). Note 
that this i?-matrix acts on the space Vj ® Vi. As shown in Eq. (IB.Sp . £j(Ai)£j(A2) is an 
operator on Vi®Vi® V2. Therefore, one can consider a product of this operator with 
-Ri2(Ai — A2). The graphical representation for the l.h.s. of the YBE is drawn as 



-Rl2(Al — \2)Ci{\i)Ci{\2 



I 

N 2 



(B.9) 



i 
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while the r.h.s. is 



A(A2)A(Ai)i?i2(Ai-A2)= 1 



(B.IO) 



[2 



Keeping in mind that blue bold dashed lines are attached to the auxiliary spaces, we shall 
henceforth draw pictures for the operators as simply as possible. Then the simplified 
graphical representation of the YBE is given by 



I 

V 2 



1 



1 



(B.ll) 



i ,2 

Using this identity repeatedly, one can show the relation Eq. f|TT]) which has the following 
graphical interpretation: 
I 



J 



1 



J 



- - - 

- - - 

1 2 



(B.12) 



Note that Eq. (fTT|) corresponds to the case where i = and j = 1. 



Appendix B.3. Recursion relations in terms of graphs 

In Sec. 3, we have derived the recursion relations for the matrices Dn and Cn- Those 
recursions are easy to understand if we use graphical representations. The recursion 
relation for Dn (see Eq. fl20l) ). 



D 



n+l 



1 

h{\n+i) 



can be expressed as 



n+ 
n _ 

1 



it 



n+ 
n _ 

1 



n+ 
n _ 

1 



it 





c(A„+i) 



n+l' 



n+r 



■X 



hi 
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Here, for simplicity, we have abbreviated the indeterminants Ai,...,A„, and A„+i. One 
can also express the recursion for C.„, Eq. f l2T]) . in a similar manner: 



n+1 
n _ 

1 



n+1 
n 



ii 



n+1 
n"' 



2 
1 



n+1 



■X 



n+1 



ii 



The matrix D„ is a lower triangular matrix while C„ is an upper triangular one. This 
can be easily seen by representing them 2x2 matrices whose elements are graphs: 



I n+lt 



D. 



n+1 



n _' 



1 



n+ 
n 



it 



1 _ 



n+_ 

_ _ n J 
^ _ iC 
i i 



it 



Cn 



+1 



1 / 



n+1 
n _ 
iC 







n+1 
n _' 

1 . iC 
Ai 

1 

iC 



[V-n + l] 



Ai 



(B.13) 



1 / 



[V-n+l] 



Appendix B.^. Mapping to five-vertex model 

As we have seen in Sec. 4, the triangular matrices Dn and C„ are transformed into T>n 
and Cn by the invertible matrix F„. This similarity transformation corresponds to a 
mapping from the six-vertex model to a five-vertex model. This relation can be easily 
understood using graphs. In the new basis, the matrix P,„ is diagonal and its recursion 
Eq. ( 146|) can be represented as 



n+ 



n 

1 _i — 







n+ 



1^ 



n 4 

1 _i_ 



v 



with Vr. 



A 

n _i_ 
i 

iJ.i 
1 i 

I 

Ai 



(B.14) 



where the vertical lines are modified to indicate the transformed operators. Since the 
recursion f H6|) uniquely determines the graphical representation for P„, all the arrows 
on the vertical line are Note that a horizontal dashed line without arrows acts on the 
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space {V) diagonally. In the main text, the matrix C„ is decomposed as C„ = ^"=1 Cn\ 
From Eq. f l5^ . we see that Cn^ flips only the state in Vj from | — t-) to | <—). Therefore, 
ciP can be written as 

f 

n _i _ _ 

1 
ii 

Since the action of Cn^ is diagonal in except for Vj, the arrows on the vertical line 
are specified as follows: 

1 

j{<-\Ci'^\-^),= i • (B.16) 
ii 

From the graphical representation, Cn^ can be regarded as a kink in J^. It is consistent 
with the picture that A^j^n becomes a Jordan- Wigner fermion in Jif when A = as 
we discussed in Sec. 4. Again we note that the the horizontal line without arrows acts 

i 

diagonally on the space. In Eqs. ( IB. 141) and ( IB. 161) . the vertex •►—--•does not appear, 

which means that the original six-vertex model becomes a five-vertex model after the 
similarity transformation. Comparing the graphs with Eqs. (!50|) and (!54l) . we find the 
weights for 6 vertices: 

A f A 



1 A 

-|-^= f{Zj,...,Zn)c{Xj), ^-j— .= 



(k) 

where j and k denote the horizontal line and the index for Cn , respectively. Note that 
the weight of the 4th vertex is determined in a non-local way, i.e., if the intersection of 

1 

the vertical line and kth horizontal line is -•-! then the weight for the intersection of 

i 

the vertical line and the jth horizontal line with j > k is given by Skj{zk, Zj) when both 
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the arrows on the horizontal edges are The weight for the 5th vertex also depends 
nonlocally on j + l,j ' + 2,...,n. In exchange for this non-locality, the combinatorial 
complexity is greatly reduced in the five-vertex model. In the original matrices Dn and 
Cn, there are 2"~^ possible configurations of arrows (t and J,) on the vertical line. After 
the similarity transformation, however, the configuration for P„ is uniquely determined 
as has seen in Eq. flB.14p . The number of possible configurations of arrows for C„ is n 
since each Cn ■* is uniquely determined (see Eq. flB.16|) ). 

Finally, we represent Z^^^'^ in terms of the transformed graphs. Recalling that 



the DWBC is not altered after the similarity transformation, i.e., Qr, 
and FnF~^ = 1, Eq. (IB. 71) is written as 



Qr. 



7DWBC 



= Tr^' [F'^QnFn (F^^C„(Ai, A„)F„) ] 
= Tr,^.[Q„Cn(Ai,...,AO"] 
= (^|C„,(Ai,...,A„,)"K). 

In short, what we have shown is graphically represented as 



(B.17) 



7DWBC 





(B.18) 



i 

with •►-!- 



I 



0. Again it should be stressed here that a merit of the five- vertex 



model introduced is the reduction of possible configurations involved in the calculation 

of ZPWBC_ 
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